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A new class of deterministic dynamical systems, termed semipredictable dynamical systems, is 
presented. The spatiotemporal evolution of these systems have both predictable and unpredictable 
traits, as found in natural complex systems. We prove a general result: The dynamics of any 
deterministic nonlinear cellular automaton (CA) with p possible dynamical states can be decomposed 
at each instant of time in a superposition of N layers involving pQ, pi,... Pn-i dynamical states 
each, where the Pken, k G [0, A" — 1] are divisors of p. If the divisors coincide with the prime 
factors of p this decomposition is unique. Conversely, we also prove that N CA working on 
symbols po, pi,--- Pn-i can be composed to create a graded CA rule with N different layers. 
We then show that, even when the full spatiotemporal evolution can be unpredictable, certain 
traits (layers) can exactly be predicted. We present explicit examples of such systems involving 
compositions of Wolfram’s 256 elementary CA and a more complex CA rule acting on a neighborhood 
of two sites and 12 symbols and whose rule table corresponds to the smallest Moufang loop Mi2(53, 2). 
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I. INTRODUCTION 


Complex natural systems exhibit a coexistence of predictable and unpredictable traits in their 
dynamical evolution. This feature seems to have been excluded from mathematical investigation, 
perhaps because of the difficulty to state the problem rigorously and because of its suspected general 
mathematical intractability mm- In this article we present a pathway to investigate the problem and 
in doing so discover an inhnite family of systems termed semipredictable dynamical systems whose 
dynamical trajectory can be sharply separated in the ordinary sum of an exactly predictable part 
and an unpredictable one, such that both parts can be neatly extracted from the trajectory at any 
time. We show how a decomposition of the dynamics can be carried out for any dynamical system 
modeled by a cellular automaton (CA) in order to extract these layers and make predictions. If the 
system is totally unpredictable, our method leads, in any case, to a reduction of the dynamics to 
independent layers that can be treated separately. These layers can then be superimposed to get the 
full solution. In many cases, however, although the full solution is unpredictable, there are nontrivial 
evolving properties of it that can be predicted at any time and we can give an explicit solution for them. 

In this article we say that a dynamical system is unpredictable if it is both deterministic and 
autoplectic in the sense given by Wolfram to this word in [3|: Autoplectic systems are spatially ex¬ 
tended systems that starting from a simple (non-random) initial condition are able to generate intrinsic 
randomness by themselves so that no regular pattern is discernible in their spatiotemporal evolution. 
Autoplectic systems include all computationally irreducible systems for which computation needs to 
be followed over time to know the outcome of the dynamical variables of the system (for such systems 
no shortcut for the trajectory generally exist). We call a deterministic spatially extended system pre¬ 
dictable if it is not autoplectic, its spatiotemporal evolution leading to discernible patterns and regular 
structures that are amenable to mathematical investigation. Most of Wolfram 256 elementary CAs 
are predictable. Known exceptions are, among others, rules 30, 110, 54 and their class equivalents 
under global complementation and left-right transformation HIM!. The complex Ginzburg-Landau 
equation with parameter values set in the Benjamin-Feir unstable regime leads to turbulence even for 
simple initial conditions and, hence, it constitutes, for such parameter values, another example 

of autoplectic deterministic (i.e. unpredictable) dynamical system. We find here however, that there 
also exist a broad class of dynamical systems of a multilayered nature, some layers being autoplectic 
and some others being predictable (even trivially predictable) and such that all layers propagate and 
interfere with each other in all locations during the spatiotemporal evolution of the system. Indeed, we 
conjecture that many CA (if not most) with a large number p of dynamical states are semipredictable. 

CAs constitute an important tool in the study of complex systems BSlElinHIS]. They have been 
satisfactorily used to model physical chemical |18j and biological m systems and often constitute 
an alternative to continuous models described by partial differential equations [20] . Although physical 
systems in nature frequently involve many microscopic degrees of freedom, a coarse graining of them 
can also be generally described by CA [211122] and the latter can thus be used to capture the collective 
behavior of pattern forming systems. Coarse grained models of CA have been shown to lead to 
predictable patterns [211 ES] even when the microscopic dynamics is unpredictable. Since coarse 
graining always implies that dynamical details are lost and these may be crucial in certain applications, 
we do not perform here any coarse graining and deal with the spatiotemporal evolution exactly. We 
show that, while keeping with all spatiotemporal information, it is still often possible to make exact 


2 


predictions of complex systems whose overall spatiotemporal evolution is unpredictable. 

Despite involving only a finite number p of dynamical states and the interactions having a finite 
range p, already the most elementary CA can already generate the highest possible complexity [IISIE]- 
In [5] we found a universal map for CA which does not depend on any freely adjustable parameter 
and which allows any deterministic CA with any number of dynamical states to be described. In 
[6] and [7] we discussed some symmetries of the universal CA map [6] getting insight in the kind of 
symmetry breaking that leads to complexity of CA evolution [7]. In this article we further explore 
these results by making use of a digit function [231 [M] that allows us to express the universal map 
for CA dynamics in a closed compact form m- One major advantage of this alternative formulation 
is then made clear, since it allows the CA dynamics to be decomposed into simpler layers that can 
be separately analyzed, the spatiotemporal evolution of the CA being a superposition of these layers. 
The outline of this article is as follows. In Section [I|] we introduce the digit function and derive some 
results in which our approach is based. Then in Section m we apply these results to CAs reducing 
their dynamics on p symbols to a superposition of (generally nonlinear) dynamical layers, each acting 
on a reduced number Pk^n of symbols, where the p^s are divisors of p. When these divisors are the 
prime factors of the canonical decomposition of p, the decomposition of the original CA rule is also 
unique. Our main mathematical result is, thus, a general decomposition property, that is contained 
in Eqs. (21) and (27) below. This result leads us to introduce the concept of p-decomposability, 


of relevance to most CA rules in rule space, and the associated concept of graded CA rules, levels 
and layer CAs. Then, we present several examples of the decomposition giving a detailed account of 
some semipredictable dynamical systems that are neither totally predictable nor totally unpredictable 
and analyze it through the mathematical methods presented in this paper. Numerical computations 
conhrm our analytical results. 

II. THE DIGIT FUNCTION 


The digit function, for p G N, A: G Z and x G M is defined as 


dp{k,x) = 


_p'^ 


^fc+l 


( 1 ) 


and gives the fe-th digit of the real number x (when it is non-negative) in a positional numeral system 
in radix p > 1. If p = 1 the digit function satishes di(A:, x) = di(0, x) = 0 and it does not relate to a 
positional numeral system. 

With the digit function we can express any real number as 


Uogp 1^1 J 

X = sign(x) p^dp{k, 


X 


( 2 ) 


k=—oo 


The digit function is a staircase of p levels taking discrete integer values between 0 and p — 1. Each 
time that x is divisible by the ascent of the staircase is broken and the level is set again to zero 
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and a new staircase begins [21] . For n and m nonnegative integers the digit function satisfies 

dp{k,x + = dp{k,x) (3) 

dp{k,p^x) = dp(0,x) (4) 

dp(0,dp(/c,x)) = dp{k,x) (5) 

dp(0, n + dp(0, m)) = dp(0,n + m) (6) 

dp(0, ndp(0, m)) = dp(0,nm) (7) 


d„p(A:,x) = dp(k,^^+pdn{k,^^^ 


= dn[k,^] + ndp ( k 


fc+i 


n 


( 8 ) 


These relationships are all easy to prove. Eq. Q an Q follow directly from the definition. Eq. 
(j^ follows by noting that dp(/c, x) is an integer 0 < dp(/c, x) < p — 1. Then 

dp{k,x) 


dp(0, dp(/c, x)) = dp{k, x) -p 


p 


= dp{k,x) 


(9) 


since [dp{k,x)/p\ = 0. Eq. @ is proved by using Euclidean division, since we can always write 
m = ap + b with a, b integers and b = dp(0, m). Therefore 

dp(0, n + dp(0, m)) = dp(0, n + ap + dp(0, m)) = dp(0, n + m) (10) 

where Eq. ^ has also been used. Eq. 0 is also proved in a similar way, using the distributive 
property of ordinary addition and multiplication as well. By using the definition, we have 


dp (k, +pdn (^k,- 


,k+l 


pknk 

X 

(np)^ 


-p 


pk+lp^k 


+ P 


X 


pk+lripk 


np 




np 


(^p)fc+l 


= dnpik,x) 


( 11 ) 


which proves Eq. Q. 

The following result is the cornerstone of our method to decompose CA dynamics as explained in 
the next section. Let po, pi,... Pn-i £ N. We have 


iV-l 


dpopi...pjv_i {k,x) — dp^ I k,- 


h=0 


h-1 

Y Pm 
m=0 


( 12 ) 


This result can easily be proved by induction. Eor = 1, Eq. ( |12| ) is trivially valid and for = 2 it 
reduces to Eq. Q. Let us assume the result valid for N factors. Then, for A/" + 1 factors we have, by 
using Eq. (fsl) 


^POPI-.PN 

N-1 

h=0 

N 


(k, x) = dpop,...pj^_^ ^/c, + pQpi... pN-idpj^ ^/c. 


{poPl . . .pAr_i)^+l 


x/p‘ 


h-l 




N-1 


k, 


(nt=oP-)^+nnnViPn) 


N 


h=0 


Pm + d. 

m=0 

h-l 

n 


k, 


X 


inYlPmr^\ 


7V-1 

n p^ 

m=0 


m=0 


which proves the validity of Eq. 
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III. THE UNIVERSAL CA MAP AND THE DECOMPOSITION OF CA DYNAM¬ 
ICS: GRADED CA RULES AND SEMIPREDICTABLE DYNAMICAL SYSTEMS 


We now establish the way in which the digit function can be used to describe CA. Let us consider 
a ID ring containing a total number of Ng sites. An input is given as initial condition in the form 
of a vector xq = (xq, ..., Each of the Xg is an integer in [0,p — 1] where superindex j G [1, W] 

specihes the position of the site on the ID ring. At each t the vector xt = (x^, specihes the 

state of the CA. Periodic boundary conditions are considered so that = xj and x° = x^ . Let 

be taken to denote the value of site j at time step t + 1. Formally, its dependence on the values 
at the previous time step is given through the mapping x^^^ = ^Rp{x{'^\ ■■■xi, ...,x{~^), which we 
abbreviate as xl_^_^ = ^Rp{xl) [5] with the understanding that the function on the r.h.s depends on all 
site values within the neighborhood, with range p = I + r + 1, which contains the site j updated at the 
next time {I and r denote the number of cells to the left and to the right of site j respectively). We 
take the convention that j increases to the left. The integer number n in base 10, which runs between 
0 and indexes all possible neighborhood values coming from the different conhgurations of 

site values. Each of these conhgurations compares to the dynamical conhguration reached by site j 
and its r and I hrst-neighbors at time t and given by 


^t= ^ (13) 

k=—r 

We will refer to this latter quantity often as the neighborhood value. The possible outputs for each 
conhguration n are also integers G [0,p — 1]. An integer number R can then be given in the decimal 
base to fully specify the rule as 

p’’+'+i-i 

R= Y, anP^. (14) 

n=0 

This is the so-called Wolfram code of the CA dynamics. We now note that, from this dehnition 


an = dp{n,R) (15) 

denotes the output of the CA when ni = n. Because n = Ylk=-r (i-®- ^ 
represented in terms of its radix p expansion), this latter equation is similar to 

«n = dp + r, ra), r\ (16) 

\k=—r / 


Thus, the universal map for CA is simply given by x^^.^ = a^j, which, by noting that x{'^^ = dp{k + 
r,nl) from Eq. (13), is given by 


X- 


U, = ^r;{xI) = dp 



k+r j+k 

Xt 


,R = 


R 


IP^ 


„k+rJ + k 


P 


R 


.P 


„k + v^3+k 


(17) 


All deterministic ID first-order-in-time CA are given by Eq. 0 for any neighborhood range and 
number of symbols. With this advantageous compact form for the universal map we do not need to give 
any table of configurations as input for the CA but just only the Wolfram code (a nonnegative integer 
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FIG. 1: Spatiotemporal evolution of the rule = da 3^+’"741831466113^ starting from an arbitrary initial 

condition on a ring of 400 sites and 300 time steps. Time flows from top to bottom. Periodic boundary conditions are used. 


Z + r+l 

R £ [0,pP ]) and the neighborhood parameters p, I and r directly. For example, the celebrated 

Wolfram’s rule 110 takes the simple form 


J _ 

^t+i ~ 


ds ( Y, no 


(18) 


\k=-l 


The maps for all Wolfram’s 256 CA with I = r = 1 and p = 2 are thus simply obtained as 

/ 1 \ 

ik+ri+k 


J _ 


dd E 2‘ 


''t ) 


R 


(19) 


Kk=-1 


by giving an input R £ [0,255]. As a further illustration, the ‘three-color’ p = 3 CA rule = 


d 


3 (EL-1 741831466113) is plotted in Fig. starting from an arbitrary initial condition. 

We can now take profit of the results derived in the previous Section for the digit function. Let the 
number of symbols p have unique prime decomposition into factors po,... ,pn-i as p = popi.. -Pn-i = 


Wh=o Ph- Then, by using Eq. (12) in Eq. (17) we obtain 
iv-i / i 


J _ 




k+r j+k 


R 


t 1 


h=0 

N-1 


\k=—r 


Z] Z p‘ 


k+r^j+k 




T 


t 1 


h=0 


, k=—r 
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N-1 


E <i» I E 




Y''- n‘ 

Z^k=-rP 


,fc + r + 


h=0 

N-1 


1 k=—r 
I 


h-1 

Pm 


E <i» I E 


(lln^O P^J ~ Y\m=oPm / m=0 

h-1 

Pm 


V' n‘ 
j^^l^k=-rP 


,fc + r + 


h=0 


: k=—r 


which, by means of Eq. Q reduces to 

iV-l 

xi+i = ^Rp{xl) = 


h=0 


pl.k=-rP t nm=0P™/m=0 


R 

riLL-rP^+^^t'^ n'^-1 


( 20 ) 


h-1 


-rj / 

m=oP^'' ^=0 


n 


( 21 ) 


We say that a CA rule ^Rp is p-decomposable if p is composite, i.e. not a prime nnmber. If a rule is 
p-decomposable, its decomposition is given by Eq. (21). Note that if p is prime Eq. (21) rednces to 
just only one term in the sum and nothing is gained. 

We observe that Eq. (21) constitutes a decomposition of the original CA on contributions coming 
from N digit functions, which represent a spatiotemporal dynamics (layer) acting on ph symbols each. 
This decomposition is thus produced by adding together all N layers, each of them being accompanied 
by a factor \^^2oPm- Such weights are responsible for the separation of the original CA into N layers. 
The latter can display complex behavior since they are generally coupled in a nonlinear fashion and the 
value on a specific layer generally depends in a complex way on the inputs from other layers. In many 
cases, however, the spatiotemporal behavior of certain layers can be exactly predicted, as we shall 
discuss in the next section through examples. First, we note that, for two given integers a > 1 and 
b, the digit function satisfies da(0, 6) = 0 if and only if a divides b. This last equation is traditionally 
expressed 5 = 0 mod a. Now, we observe that, indeed, from Eq. (21) (since 0 < dp^(0,y) < Ph — ^ 
regardless of the value of the quantity y) 


dpo(0,x^+J = dpo(0, 'i?;;(x^)) =dpo ( 0, 


R 


w' „k+r^3+k 
p2^k=-rP 


dpopi (0, x^_l_i) — dpupj(0, Rp{xi)) — dpQ I 0, 


R 


P- 


+ podpi 0, 




,k + r^j + ><^ 

R 


Y -pA^k=—rr Pq 

‘^POPlP2(0) ^t+l) = ^P0PlP2{^^ Rpi^i)) ~ WE k + r j + k J 

\ pl^k = -rP J 



W* „fc+rE + '= 
pA^k^-rt' t PqPi^ 


etc. 


And this system of equations, no matter how many layers we have in the composition, can always be 
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solved for the layers as 

■ ^0, - 


■*P0 I 

V p- 


,fc + r„J + '= 


d 


■Pi 


0 , 


d 


■P2 


0, 


P 


Efc=_rP''+’'a:-t 

R 

p2^k=-rP t 

R 


= dp,io, ^R;ixl)) 

dpopi(0, ^K^pixD) - dp^{0, ^Rlpixl)) 


= d 


Po 


■Pi 


0 , 


^R^ixl) 


dpoPiP2(0> ^p(a:t)) -dpopi(0, ^RlixD) 

- - - - - = d 


PoPi. 


PoPi 


■P2 


0 , 


Po J 
PoPl 


. . etc. 


( 22 ) 


where Eq. ([^ has been used in getting to the last identities. We then observe that Eq. (30) can be 
written as 

^RUxi) 


N-l 


^t+1 


= ^r;{xI) = dpj 0, 


n h—1 

m=0 Pm ^ 


h=0 


h-1 

n Pm 

m=0 


(23) 


which is an alternative form of Eq. (21) and which, indeed, follows directly from Eq. (12). Thus we 
have the N equations given by 


0 , 


X' 


t+i 


Ph I Trh-1 

V 1 lm=0 Pm , 


= d 


0 , 


R 


I ’ yi p'^+^xi+’‘ T-ih-i 

\ P ‘ nm=oPiii. 


/l = 0,1, 


,iV-l 


(24) 


which govern the spatiotemporal evolution of each layer h in the composition as a function of the 
p-decomposable rule ^R!p{x{) with p = O^o^P^- that, in general, as a consequence of 

Eq. pll) 




= dpo +Podpi 0, — TpoPid 


x-l 


N-l 


Po 


■*P2 


POPl 


+... - dp^ 


0, 


XI 


h-l 


n h—1 j XX 

m=0 P'^ / m=0 


/i=0 


\{Pm (25) 


at every j and t. Therefore, if we dehne yf= dp^ ( 0, „h-i -) so that takes values only on the 

_ \ llm=oP"* / 


subset [0,p/i — 1], Eq. (25) can be rewritten as 

3 Oj I Ij I 2,7 I \ ^ h,j I I 

Xt = Vt +Poyt TPOPlPi + • • • = 2^ Pt 11 Pm 


N-l 


h-l 


h=0 


m=0 


and, therefore, Eqs. (24) take the form 


h,j _ 


Vt+i 


= d 


Ph 


0, 


R 


■wi „fc+r'wJV-l h.j + fcprTi-l -nh-l / 

p^k^-rP yh^oVt nmJoPm/ 

If we compare this latter expression with Eq. (by using also Eq. (Q there) 

R \ 


h = 0,l, 


^t+i 


= do 0 




p 


(26) 


(27) 


(28) 


we observe that the layer h is described by a CA which is independent of the other layers only if the 
following relationship is satished 


_ 


Pt+i 


= d 


Ph 


Or 


P 


R 

nt\\Pm Tih-i 


= d 


n n—i 

m=0 Pm / 


■Ph 


0 , 


lW 


E l ^k + r h,j + k 

k=-rPh Vt 

Ph 


(29) 
























for some integer L^^'> such that 0 < — 1 and for all possible values of G [0,Ph ~ !]• 

If a similar expression is satisfied for each layer h, we say that the rule ^Rp is a graded rule. In a 
graded rule all layers are decoupled and described by a CA. Graded rules are thus a subset of the 
more general p-decomposable rules. 

Eqs. (21) and [21) are the main result of this article. To understand the usefulness of these 
expressions and to better grasp the concept of a graded CA rule it is appropriate to think as well in 
the opposite process of constructing a p-decomposable rule out of the composition of more elementary 
CA rules. The following text and the examples in the next section should be helpful to comprehend 
this abstract approach and its associated unfamiliar notation. 

Let us construct a graded rule out of more elementary CA. Let thus p 2 ) ^tc. 

be CA rules acting on alphabets of po, pi, p 2 , ... symbols each. Our above approach means that we 
can construct a p-decomposable CA rule ^Rp, with p = poPiP 2 ■ ■ ■ out of these CA rules in the following 
way 


4+1 


= ‘R'.ix’) = +4; (d„ (o,rc;))+j.„ 


+. .. 



where {k G [—r,/]) and xij^i are all G [0,p — 1]. Each denotes an elementary CA rule 

acting on p^ symbols. All rules constructed through Eq. (30) are graded CA rules. In Eq. (30) we 
say that rule is the ground rule, that rule is lifted by po and, in general, that rule 


is lifted by a factor pQpi .. .ph-i = Y[m=oPm- general, we call po, popi, P 0 P 1 P 2 , ■ ■ ■ the levels and 
the CA rule on each level the corresponding layer CA. Thus, we note that, from Eq. (p0|, Eqs. 

(27) take the form 


0 , 


X 


t+i 


= ^Tfh) r 
Ph 


0, 


Ph I Y~fh—1 ^ I Ph \ '~^Ph \ y-fh—l ^ , , 

\ 1 lm=0 Pm./ \ \ 1 lm=0 Pm J J 

and, we observe that, for every j and t, this latter equation becomes 


h = 0,1, 


,iV- I 


( 31 ) 


*1 = ‘i'T (nl’ 


/i = 0,1, 


,A^-1 


(32) 

which is simply the equation of a decoupled CA acting on p^ symbols. The y^’^A evolve within 

the graded rule and are easily extracted from it through their definition y’f'^ = dp^ [ 0, „h-i - 

\ llm=0P"i 

Conversely, the dynamical state x{ of the graded rule is constructed out of the yf’'^ ’s through Eq. 

(^. 

Each layer CA ^L^^p/ has Wolfram code where l^} = dp{nh,L)^)) is an 

integer G [0,p/i — 1]. Now, by looking at the derivation of Eq. 0 from Eqs. ( [T^ to ( [l6| ), we observe 
that, from the construction, the graded CA rule ^Rp has Wolfram code 


„l + r+l 


-1 


„l + r+l 


-1 


N-l 


h-l 


^Ph 


R= Yi ?"“»= E f>"E'i: 

n=0 n=0 h=0 

and, thus, its coefficients in the rule table are 


Epf 


fc+r 


k=—r 



m=0 


N-l 


I 


a^f 


— dp(n, R) — ^ dp^ ( ^ pjf^ dp^ ( 0, 


h=0 


\k=—r 


dp^ (fc + r, n) 
l\m=eiPm , 


h-l 


UPrn 


(33) 


(34) 


m=0 
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for n = 0,... — 1. 

The following is a summary of some observations: 

• Out of arbitrary CA rules ^ 


, we can construct a CA rule 


^Rp which is a graded rule (superposition) of the former. This is provided by Eq. (30). 

From the knowledge of the graded rule = ^Rpop^p 2 ... alone, the layer CA rules , ^L^^pi 


/r(2) r 
^ P 2 > 


can be extracted at every time and every location. This is what Eqs. (31), resp. 


Eq. (32), mean. 


• If all factors of a graded rule are equal, po = ... = pn-i then, although a permutation of the 
layer CA rules do not change the layer content of a graded rule, the spatiotemporal evolution 
of the resulting graded rule is distinct to the original one (before the permutation) for the same 
initial condition. Thus, A^! same-layer-content graded CA rules with different spatiotemporal 
evolutions (when starting all from the same initial condition) can be constructed (the layer CAs 
being the same in all cases, but acting on different levels of the corresponding graded CA rule). 
The levels break the permutation symmetry of the layer CAs within the graded rule. Indeed, a 
graded CA rule with levels I, po, p ^,..., Pq~^ has the form of a Galois resolvent (see [5^, pp. 
35-37 for an introduction to this concept and to its usefulness in group theory, as well as Galois’ 
memoir in Appendix I of [25j, i.e. the function V on Lemma II). 


• We have seen that any factors po,. .., Pn-i give rise to a decomposition in levels I, pO) PoPi; • • •) 
Y\h=o Ph- Thus, if some factors are repeated, let aj denote the multiplicity of the factor pj and 
let M denote the number of distinct factors. Let the index j run over these distinct factors. We 
then find that there are 


iV! 


n M — l I 
7=0 


(35) 


different ways of arranging the same factors into levels. This again leads to respective distinct 
decompositions with same-layer CA content. For example, for p = 30 and factors (2,3,5) there 
are 6! graded rules (since we have ao = ai = 02 = 1 and hence M = N = 3) that can be derived 
by arranging the factors in different levels (I,Po,PoPi) as (1,2,6), (1,2,10), (1,3,6), (1,3,15), 
(1,5,10), (1,5,15). To these levels correspond, respectively, the arrangements of the layer CAs 
with (po,pi,P2) symbols given by (2,3,5), (2,5,3), (3,2,5), (3,5,2), (5,2,3) and (5,3,2). 


Eq. 


is not equivalent to Eq. (21) but a specific case of the latter. We have that a graded 


rule is always p-decomposable but the contrary is not true. Indeed, if p is composite there are 

I I T' I 1 

pP rules that are p-decomposable of which only a subset are graded rules: A p-decomposable 
rule is a graded CA rule only when each layer’ dynamics, as generally given by Eq. (27), can 
be brought to the simple form of a layer CA given by Eq. (32). 


• Although there is a unique decomposition of a p-decomposable rule into the prime factors 
of p, there are a number of other possible decompositions in levels which equals the total 
number of multiplicative partitions [26] of p. The method above is valid to construct any such 
decompositions. For example, since p = 30 has five multiplicative partitions 2x3x5 = 2x15 = 
6x5 = 3xI0 = 30a CA rule with p = 30 admits five possible decompositions in factors 
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(2,3,5), (2,15), (5,6), (3,10), (30), each describable by Eqs. (21) and (27). If all N factors of 
p are distinct (i.e. if p is squarefree) the number of possible decompositions of the graded CA 
rule is given by the Bell number B]\f, which can explicitly be calculated as m 


N 


Bn = Y1 


k=0 


N 


1 


k=0 j=0 


(36) 


where | ^ | is a Stirling number of the second kind. 


Compared to a general p-decomposable rule, the layer dynamics greatly simplifies in the case of a 
graded CA rule, since each layer is independent of the others (even when all interfere with each other 
in the spatiotemporal evolution of the full graded CA rule). Although the layers are all entangled in 
the graded CA rule, they can be separated with help of the digit function, as established above. 

The above conclusions have still the following important implication. If any layer in Eq. (21), 
displays a complex and unpredictable behavior the compound CA rule ^Rp will also generally display 
a complex and unpredictable behavior. However, certain layers within a p-decomposable rule may be 
predictable, even when some others are not (and the composition is certainly not). We term such 
CA rules semipredictable to emphasize the fact that shortcuts or even closed forms can sometimes be 
found for the orbits of some of the separate layers. The examples in the next section substantiate this 
statement. 


IV. EXAMPLES AND DISCUSSION 


A. Graded CA rules from Wolfram’s elementary CA 


Let us first consider an example of a construction of a graded CA rule out of two simpler rules. 
Wolfram’s rule ^30^ is well-known to display an unpredictable, chaotic behavior [28]. This rule belongs 
to Wolfram’s 256 elementary CA with two symbols (p = 2) and a neighborhood with one left and one 
right neighbors (Z = r = 1). On the contrary, Wolfram’s rule ^1502 is a predictable (albeit complex and 
subtle) |29| CA displaying a nested pattern. Wolfram’s rule ^1502 has map +2 x{ +2 x{~^ 

where +2 denotes addition modulo 2 (ISllT]) and, thus, is a linear CA. Such CA are not autoplectic 
|3|. The explicit solution for the orbit as a function of the initial condition can indeed be found in 
some references (see e.g. m)- Although such solution depends on a sum with a number t of terms 
taken from the initial condition and weighted with binomial coefficients, fruitful use can be made in 
most cases of the Lucas’ correspondence theorem [7|, to simplify the resulting expression. We shall 
discuss linear CA in more detail elsewhere in connection with Lucas’ theorem. 


We can use Eq. (30) to construct a CA with p = A, pQ = pi = 2 out of these two rules as 


J _ 

^t+i ~ 


^Rl{xl)= ^302 (^d2 -h2 • ^ 502 ^ 2 ( 0 , 


d2 ( ^ 2^+id2 (o,x^+^) ,30 I +2d2 ( ^ 2^+M2 ( 0, 



^fe=-i 




,150 


(37) 


^fc=-i 
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This is a graded CA rule acting on four symbols 0, 1, 2, 3 which are the values that can take at 
each time t and for every j. The Wolfram code of the graded rule is easily calculated from Eq. (33) as 

63 


43-1 

R = = 


(38) 


n=0 


n=0 


where the a^’s are given by Eq. (34) as 


— 


d 4 (n. R)=Y: { E 2*«d. (o. ( 

h=0 \k=-l ^ ^ / 


n = 0, ...,63) (39) 


with = 30 and = 150. Thus, the rule vector of the graded CA rule is obtained from this 
expression 

(ao, ai,..., oes) = 

(0, 1, 2, 3, 1, 1, 3, 3, 2, 3, 0, 1, 3, 3, 1, 1, 1, 0, 3, 2, 0, 0, 2, 2, 3, 2, 1, 0, 2, 2, 

0, 0, 2, 3, 0, 1, 3, 3, 1, 1, 0, 1, 2, 3, 1, 1, 3, 3, 3, 2, 1, 0, 2, 2, 0, 0, 1, 0, 3, 2, 

0, 0, 2, 2) (40) 

The rule has a huge Wolfram code of order R ~ 4®^. The spatiotemporal evolution of the rule and 
the layer CAs are dictated by Eqs. (32) and (26), which take the form 


2/m 

= '30i( 

= d2 ( 


0,1 

2/0 = 

= d2(0. 

Xq) 

(41) 




\k=-l J 

( ^ \ 





2/m 

= H50^ 

= d2 

^ 2^+’'2/,^■+^150 

1,1 

2/0 = 

= d2(0. 

4/2) 

(42) 




\fc=-i / 





j 

Xt+l 

o,i , 

= 2/t+i + 

22/iYi 





(43) 


starting from an initial condition Xq. Asa consequence of the main result of this article, these equations 
yield the full description of the spatiotemporal dynamics of the graded CA rule with 4 symbols given 
by 


= ^Rl(xi) = dd 






(44) 


^fe=-i 


with R given by Eqs. (38) and (40). 

We have performed numerical calculations of the spatiotemporal evolution of the CA above. On 
one hand, it is straightforward to obtain the evolution of the rule from Eqs. (41) to (43). On the other 
hand, calculating the evolution directly from Eq. (44) is a bit more demanding because of the large 
powers of the radix p = 4 in the digit function. This problem is computationally overcome in the light 
of our previous work |5]: We consider the rule vector, as given by Eq. (40), and use the universal map 
for CA derived in [5] 


pr+l + l_i 


x^+1 = ^R^pixl) = anBln- 


^k+r j+k 

P , 


(45) 


n=0 


k=—r 
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FIG. 2: Spatiotemporal evolution of (a) (b) = d 2 ( 0 ,a;|) and (c) = d2(0,3^j/2) all calculated by iterating Eq. ( |44[ | 

(see text), for the CA rule with code given by Eqs. ( |38[ l and ( |40[ l and for initial conditions corresponding to = 3 and = 0 
for j 7 ^ 250. We observe that the rule is four-color graded rule, with layer CAs consisting of the two-color Wolfram rules ^30^ as 
ground rule, and ^1502, as shown respectively in (b) and (c). Time flows from top to bottom in each panel. Periodic boundary 
conditions are used, j E [1,500] and time t E [0,240]. 


where 




1 ( x + y 


x-y 


= 2 (sign(x + y) - sign(x - y)) 


(46) 


\x + y\ \x-y\^ 

is the ;B-function for any real numbers x, y. This function returns sign(y) if —y < x < y, sign(y)/2 
if X = y and zero otherwise. Eq. (45) is equivalent to Eq. 0 and overcomes the problem of large 
powers of p. While the digit function has mathematical advantages compared to the .B-function, the 
latter is more suitable for intensive _computations. In Eq. 

“ ' (b) 


the Qn given by Eq. (40). In Fig. 


(45) we thus take p = A, I = l,r = 1 and 
we show the spatiotemporal evolution of (a) x^ 


'^t ) 


d2(0, X- 


(c) d2(0,xi/2) calculated from Eq. ( |44[ ) (i.e. Eq. (45) with the rule parameters given by Eq. 
(|40[)) and for a simple initial condition. We thus see that such a CA rule allows both chaotic and 
coherent signals to be spatially transmitted together in different layers, starting from a simple initial 
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FIG. 3: Spatiotemporal evolution of (top), = d 2 ( 0 , (bottom left) an d = d 2 ( 0 ,x{/ 2) (bottom right) all calculated 
by iterating Eq. ( |44[ | (see text), for the CA rule with code given by Eqs. ( |38[ | and l |47[ l and for a random initial condition. The rule 
is a four-color graded CA rule, with layer CAs consisting of the two-color Wolfram rules ^204^ (the identity) as ground rule and 
^1372 (capable of universal computation PHIS]) as lifted rule. The graded CA rule in the top panel is thus a semipredictable 
system. Time flows from top to bottom in each panel. Periodic boundary conditions are used, j G [1,400] and time t G [0,200]. 


condition within the same compound signal. Although both layers interfere in forming the pattern of 
the graded rule, they can be neatly separated by means of the digit function into two different signals 
both employing two symbols and behaving quite differently: One of them, albeit complex, has a large 
degree of coherence; the other one is incoherent. All these statements are true for any arbitrary initial 
condition. We have calculated the graded rule both from Eqs. (43) to (43), and from Eq. ( |44| ) (resp. 
Eq. (45)), computationally confirming the equivalence of both approaches. 

Another example is provided by layer CA rules ^204^ and ^137^ with Wolfram codes = 204 
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and = 137. The rule vector of the graded rule is calculated from Eq. (39) as 

(ao, ui, • • ■, oes) = 

(2, 2, 0, 0, 3, 3, 1, 1, 0, 0, 2, 2, 1, 1, 3, 3, 2, 2, 0, 0, 3, 3, 1, 1, 0, 0, 2, 2, 1, 1, 

1, 3, 3, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 2, 2, 1, 1, 3, 3, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 2, 

2, 1, 1, 3, 3) (47) 

The spatiotemporal evolution of the graded rule is shown in Fig. [^together with the ones of its layer 
CAs. CA rule ^204^ is trivially predictable since it is equal to the identity and one has for 

every t. On the contrary, rule ^137^ displays complex and unpredictable behavior (the rule is just the 
global complement of rule ^110^ which is known to emulate a universal Turing machine mm) and no 
shortcut for the orbit of yl’^ is possible. Consequently, no shortcut is possible for the graded CA rule, 
whose evolution is given by xj = y^'^ + 2yl’^ (top panel). The unpredictable character of the graded 
rule is entirely given by that of the rule 


1137^ 


Thus the rule vector in Eq. (47) yields a graded rule 


^R\ which is a semipredictable system one layer being trivially predictable and the other one not. 

We thus note that all 4^^ = 4^^® ~ 3 • 10^® CA rules ^R\ are p-decomposable. Of these, only a 
subset consisting of 2^^2^^ = 256^ = 65536 are graded rules, obtained from choosing any two of the 
256 elementary Wolfram CA rules as ground and lifted layer CAs. 

We now construct more complex graded CA rules with p = 16, arising as the result of the compo¬ 
sition of the CA rules ^1242, ^137^ and ^1932. All these CAs are related by global complemen¬ 

tation and/or left-right transformation and all four are capable of universal computation. Thus, the 
corresponding 4! = 24 graded CA rules that can be formed with these four layer CAs are all universal 
as well. Let us consider the graded CA rule given by 


Oj 

Vtii 


Vtii 


2,1 

Vtii 


3,1 

Vt+l 


= AlO, 


424^ 




1,1 


^137^ 


493^ 


Vt 


Vt 


Vt 


3j 


= d 


= d 


= d 


= d 


The above system of equations completely specify the spatiotemporal evolution of the rule. Let us 
call it rule A. Now, if we e.g. make the transformation /i —?■ 3 — /i in y^’^, Vj and t, the result is 
another graded CA rule with the same layer content than the original one, but leading to a different 
spatiotemporal evolution for the same initial condition. This is so because of the symmetry breaking 
introduced by the different level weights 1, 2, 4 and 8. Let us call this rule so obtained rule B. In Fig|^ 
the spatiotemporal evolutions of rule A (top) and rule B (bottom) are plotted for an initial condition 
given by Xq = di6(l, j). Although the decomposition in layers of both rules is the same, as is the 
initial condition, the spatiotemporal evolution looks very different. Rule A displays a proliferation of 
coherent structures creating a kind of granular texture. Rule B seems to lead to more disorder. 





k+r„p,j+k 

yt ) 


‘yk+r l,j+k 

^ yt 5 


ryk+r 2, j+k 

^ yt ) 


yk+r ^,j+k 

^ yt ) 


no 


124 


137 


193 


Vo^ = d2(0,x/,) 


= d2(0,x/|/2) 


= d2(0,x/|/4) 


= d2(0,x/|/8) 


(48) 


(49) 


(50) 


(51) 


(52) 
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FIG. 4: Spatiotemporal evolution of for rule A (see text) given by Eqs. \ to j52| (top), and for rule B (bottom) obtained 
by a permutation of the layers of rule A (see the text). Both rules have number of symbols p = 16 and neighborhood range 
/ + r + 1 = 3. The initial condition is equal in both cases and taken as Xq = di6(l, j). Thus, although the layer content of both 
graded rules is the same and the initial condition is also identical in both cases their spatiotemporal evolutions are quite different 
because graded rules are not invariant under permutation of the layers. Time flows from top to bottom in each panel. Periodic 
boundary conditions are used, with j E [1, 528] and time t E [0, 250]. 
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B. A non-graded p-decomposable rule with semipredictable dynamics 


In the previous examples we have considered the construction of p-decomposable graded rules from 
layer CA rules. In the following example we consider the opposite process: Given a p-decomposable 
CA rule, we derive and analyze its layers and attempt to reach conclusions on its dynamics. We 
consider a p-decomposable rule that, as we shall see is not a graded rule. 

Let us consider a rule with p = 12, / = 0, r = 1, range p = l + r + l = 2, Wolfram code 

143 

= (53) 

and rule vector 


(ao, ai,..., ai 43 ) = (0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 1, 0, 4, 5, 2, 3, 7, 6, 10, 
11, 8, 9, 2, 5, 0, 4, 3, 1, 8, 11, 6, 10, 9, 7, 3, 4, 5, 0, 1, 2, 9, 10, 11, 6, 7, 8, 4, 

3, 1, 2, 5, 0, 10, 9, 7, 8, 11, 6, 5, 2, 3, 1, 0, 4, 11, 8, 9, 7, 6, 10, 6, 7, 8, 9, 11, 

10, 0, 1, 2, 3, 5, 4, 7, 6, 10, 11, 9, 8, 1, 0, 5, 4, 2, 3, 8, 11, 6, 10, 7, 9, 2, 4, 0, 

5, 3, 1, 9, 10, 11, 6, 8, 7, 3, 5, 4, 0, 1, 2, 10, 9, 7, 8, 6, 11, 4, 2, 3, 1, 0, 5, 11, 

8, 9, 7, 10, 6, 5, 3, 1, 2, 4, 0) 


(54) 


This rule provides the map xi_^^ = ^R\ 2 {xi) = x{ * x{~^ which is described by Eq. (17) and which 
shows that the CA rule ^R \2 above can be understood as a binary operator * (law of composition) 
S' X S' —>• 5 where S is the set of the integers G [0,11]. We note that, n = di 2 (0,n) + 12di2 (l,n) 
(n G [0,144]) and that n{ = 12x{ + where xi = di 2 (l,n^) and = di2{0,n{). Thus, since 
x{_^_i = a^j = xi * xi~^ = di 2 ^l,nj^ * di 2 ^0,n^^, it proves useful to arrange the rule vector on a 
table so that the output = di 2 (1, n) * di 2 (0, n) is obtained from reading the table for any di 2 (1, n) 
listed in the rows and and di 2 (0, n) listed in the columns of the table. The latter has the form 
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7 
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2 

11 
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2 

3 

1 
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5 

6 

5 

3 

1 

2 

4 
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This table corresponds to the Cayley table of the smallest Moufang loop Mi 2 (S 3 , 2) [3T] of order 12 
(see [32] for a classification of small Moufang loops). A law of composition *,SxS—)-Sisa Moufang 
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loop if it satisfies the following identities, called the Moufang identities (3^ 


z * {x * {z * y)) 

= {{z * x) * z) * y 

(55) 

X * {z * {y * z)) 

= {{x * z) * y) * z 

(56) 

{z * x) * {y * z) 

= (z * {x * y)) * z 

(57) 

{z * x) * (y * z) 

= z * ((x * y) * z) 

(58) 


for x,y,z G S. It can be checked from the table above that this is indeed the case. Here S is the set of 
integers G [0,11]. We see from the table above and the Moufang identities that the operator * is both 
non-commutative and non-associative. However, x = z, x = y or y = z, the resulting subalgebra is 
associative. We thus expect that (as a consequence of losing the associative property) a high degree 
of disorder/unpredictability in the general spatiotemporal dynamics and (as a consequence of the 
existence of associative subalgebras) the existence of some regimes where the dynamics is simpler to 
describe. 

Let us investigate the above facts with the mathematical methods presented in the previous sec¬ 
tion. We have that p = 12 = 2^ • 3 is composite and admits multiplicative partitions with factors 
{P 0 ,P 1 tP 2 ,P 3 ) given by (2,2,3), (4,3), (6,2) and the trivial (non-decomposed) one (12). We study the 
partition (6, 2) with the layer po = 6 as ground layer and the layer with pi = 2 lifted by 6. From 
Eq.(26) and Eqs. (27) we have that the dynamics of the layers and of the p-decomposable rule are 
given by 


Oj 

Vt+i 


vl+i 


a 12*^ 

Z^n=0 


— de 0, 


= da 0, 


12EL-112 ''+! eLo %n:rio p- ^ 

a 12"^ 
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6. 12EL-i12'=+ieLo% 
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n 

1 In 


1 Pm 


Vo^ = d6(0,x/|) 


2 / 0 ’^ = d2(0,x/|/6) 




= yt7i + %+i 


for the On’s given by Eq. 


After a little algebra, these equations reduce to 


143 


0,1 

Vt+i 


1,1 

Vt+i 


= dfi 


0 ,?- 
n-Vt 




= da 


0,^anl2' 

^ n=0 / 

^ T 143 

^ n=0 / 


Vo^ = d6(0,x/,) 


yl’^ = d2(0,x^/6) 


J _ 

^t+i ~ 




(59) 

(60) 
(61) 

(62) 

(63) 

(64) 


If we put nt = y^'^ ^ -I- 6yl'^ ^ -|- I2y^’^ + 72yl'^ we note that if a 


•nt 


< 6, yl^ii = 0 and if ant ^ 6 
then yl^i = 1. By looking at the above table giving the 144 entries of the Moufang loop, we see that 
yl^^i = 0 when nt takes values on the 6x6 square blocks on the diagonal of the table. Conversely, we 
have that yl^i = 1 when nt falls in any of the nondiagonal square blocks. This in turn means that 
yl^^i = 1 when either y]’^ ^ or yl’^ are equal to one, but not both, regardless of the values of y^’^ ^ 


and y^'K Thus the output yt'^i is the XOR function (addition modulo 2+2 [Si [7]) of y/ 
and we find that 


1,1-1 1,1 

or y^ 


ylti 


( 1 143 

= ddo 

V n=0 




= yY ^2yY-^ ^A^%yY ^y]:l-^) 


(65) 
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where the binary operator +2 has table 


+2 

0 

1 

0 

0 

1 

1 

1 

0 


So, the lifted layer corresponds to the simple CA with code *^62 [Z| and rule vector (ao,...,a 3 ) = 
(0,1,1,0), which is a layer CA, independent of the ground layer. This rule is predictable. For, let us 
define the shift operator T as 




( 66 ) 


it can be easily proved that the set of configurations of the CA forms a ring under addition modulo 2 
and the action of the shift operator. Therefore, since 


y\t, = +2 = (1 +2 T)y^^' 

starting from an arbitrary initial condition y^’^ we obtain the solution as 


2/t= (1 +2 T)Vo’^ = 

k=0 


mod 2 = 


E 

k=0 


Vo 


Ij+fc i 2 


= d2 0 ,^ 


k=0 


^ \ 1,1+fc—i 

kh 


(67) 


( 68 ) 

(69) 


where the binomial theorem has been used. Thus, the solution at any later time t is explicitly known 
from this equation for any arbitrary initial condition. 


The meaning of the layer described by Eq. (63) which, as we have seen, reduces to Eq. (65, with 
solution given by Eq. (69) is now clear to explain. If one considers the sets of values Sq = {0,1, 2,3,4, 5} 
and Si = {6,7,8,9,10,11} then Eq. (65) specifies at each time to which set belongs from the 
information of those of xi and xi~^. Equation (65) establishes that, at time t + 1 the value 
belongs to set Si if either xj or xi~^ belong to Si, but not both. A most remarkable fact is that we 
can exactly predict by means of Eq. (69) at every time and for every initial condition Xg, to which set 
does xj. belong. 

We find, however that the ground layer is not independent of the lifted layer and, hence, the 
Moufang rule is not a graded rule. The behavior of Eq. (62) is far more complex and no shortcut for 
the trajectory seems possible. Thus, the full evolution dictated by Eq. (61) is unpredictable as well. 
The conclusion of all this is that while we cannot predict the actual value of xj in this case, we can 
exactly predict at every time, from Eq. (69), to which set Sq or Si belongs the value of x^. 

In Fig. [^we show the spatiotemporal evolution of x{ (top), y^’^ = d6(0,x^) (center) and yl’^ = 
d 2 (l,Xj/ 6 ) (bottom) calculated from Eqs. (62) to (64) for an initial condition Xq = di 2 (l,j). We 
observe that, although for x{ (top) and the ground layer y^’^ (center) the spatiotemporal evolution is 
unpredictable, the one of the lifted layer yl’^ (bottom) is exactly predictable and given by Eq. (69) 
(we have numerically checked the correctness of this statement). This is so, of course, regardless of the 
initial condition. The contribution of the lifted layer to the overall spatiotemporal dynamics is clearly 
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FIG. 5: Spatiotemporal evolution of (top), = ( 16 ( 0 , 3 ::^) (center), and = d 2 ( 1 , 3 ;j / 6 ) (bottom) calculated from Eqs. 
| |62[ | to ( |64| l for the CA rule with p = 12, / = 0 and r = 1 with rule vector given by Eq. l |54| | and for an initial condition given 
by Xq = di 2 (l,j'). While the spatiotemporal evolution of the ground rule (center) is unpredictable, the one of the lifted layer 
corresponds to the simple CA ^62 which is predictable for any arbitrary initial condition. Time flows from top to bottom in each 
panel. Periodic boundary conditions are used, j G [1,400] and time t G [0,200]. 
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discernible. Although the lifted layer is predictable and independent of the ground layer, the latter is 
dependent of the lifted layer and thus, the p-decomposable rule governing the overall dynamics is not a 
graded CA rule. Let us try to better understand the dynamics of the ground layer, which is obviously 
correlated to the lifted layer, as Fig. [^(center) and (bottom) show. Indeed, from Eq. (62) there is 
an explicit dependence on the lifted layer through and We have now four possibilities: 1) 


,bf-i 


= 1. In cases 


Vt -Vt = 0; 2) = 0 = 1 - ; 3) = 1 = 1 - and y^'^ = y^ 

1) and 2) the action of Eq. (62) at time t reduces to a CA °Ag, in case 3) to a CA and in case 
4) to another CA The rule tables of these CAs, which act on 6 symbols are given in the tables 
below 
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Thus, we reach the following conclusion: While the lifted layer is independent, the ground layer can 
be understood as being unfolded by the action of three different layer CAs of 6 symbols that act 
locally at j and t. Which CA is chosen to act on that layer, depends on the values of yl’^ and yl’^~^ 
of the lifted layer at j and t. This is interesting because, in their course of action, complex systems 
frequently chose within alternatives as a function of other independent factors, and this process can 
be modeled by the CA above. 


V. CONCLUSIONS AND PERSPECTIVES 

In this article we have presented a method to decompose any nonlinear CA acting on an alphabet of 
p symbols into a superposition of N CAs, acting each on a number of symbols ph, h = 0,..., N — 1 that 
is a divisor of p. Since the prime factor decomposition of a number is unique (by the celebrated Euclid 
theorem) so is also that particular CA decomposition unique for each CA in rule space. This general 
result applies to all CA in rule space and has lead us to discover a new class of dynamical systems 
that we have termed semipredictable dynamical systems: In these, although the overall spatiotemporal 
dynamics is unpredictable, certain traits of the evolution can exactly be predicted. We have given 
explicit examples of semipredictable systems, from those involving elementary CA, to a rule of range 
p = 2 acting on y = 12 symbols with the structure of the smallest Moufang loop. We have seen that, 
in spite of the spatiotemporal dynamics being chaotic and unpredictable, the CA can be decomposed 
into two layers, one of them unpredictable and the other displaying a predictable spatiotemporal 
evolution/pattern for any initial condition. This is a remarkable fact since it shows that a chaotic 
scalar signal can carry with itself coherent information that can be extracted from a digital analysis 
of the signal. 

The semipredictable dynamical systems presented here may find applications in the modeling of 
biological and complex systems. Most complex systems found in nature possess both predictable and 
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unpredictable traits. For example, a certain disease on a complex organism can have a predictable 
evolution, with the dynamics of the complex organism that acts as host being unpredictable in many 
details that are relevant to the overall dynamics. By means of our model we have seen that, even 
when predictable and unpredictable layers describing the dynamics are correlated, and hence not 
totally independent, they have sometimes robust features (closed substructures) that warrant their 
separate tractability and solvability. We give a summary of certain applications that seem naturally 
suited to the mathematical approach presented here. 

• Physics: Together with the p\n fractal decomposition that we have very recently introduced [23] , 
graded CA might be applied to construct useful tensor products of wave functions in quantum 
mechanics. In fact, the connection of the digit function to Hilbert space is natural from our 
previous work |23| and from the principle of least radix economy |23| . We have recently proposed 
the latter as a possible pathway to derive physical laws in a unified manner (encompassing both 
classical and quantum physics). The idea of semipredictable dynamical systems might mediate 
the quantum-classical correspondence in an analogous manner to the behavior found in this 
article for the ground and the lifted layers of the Moufang loop example. The key physical 
quantity to be addressed is the Lagrangian action S [23| which is also connected to entropy 
IsaiMlESI and which yields a natural physical radix p (which can be thought as related to p in 
this paper) that establishes the character of the physical laws [23]. Although we suggest that 
CA methods and fractality are of major importance in quantum mechanics, we do not assume 
any discrete or fractal spacetime at the Planck scale in and we do not need such assumption 
in view also of the results reported in where we show that fractals can be constructed out 
of continuous and differentiable functions that act as ‘mother functions’. Such latter complex¬ 
valued functions can be made Lorentz covariant and so shall then be any fractal objects derived 
from them through the method presented in |24j as well. Other approaches based on CA 
[36H4U| may also benefit from the mathematical methods introduced in this paper. Another 
application might be found in the analysis of complex fluids by means of lattice Boltzmann 
models [a [IT]: Lattice fluids involving many dynamical states can be decomposed in layers 
with the methods presented in this paper to get insight in their interactions and how the latter 
affect their spatiotemporal evolution. 

• Mathematics: By considering cellular automata of neighborhood size p = 2, the abstract al¬ 
gebraic properties of closed binary operators, also called laws of composition or magmas m, 
can be investigated by studying the spatiotemporal evolution of their associated CA rules. The 
Moufang loop presented here is an example. Finite groups and any other magmas [HEMS] 
can also be investigated in a similar manner. The idea of p-decomposability of the CA dynamics 
presented here is directly related to the Lagrange theorem of finite group theory. CA layers in 
graded CA rules are analogous to normal subgroups of finite groups. A graded CA rule can also 
be understood as the dynamical counterpart of a direct sum of discrete operators. Of particular 
interest is also the connection established with multiplicative partitions |26) . Thus, our formal¬ 
ism provides an alternative pathway for the computational investigation of abstract algebraic 
structures. 

• Complexity science: Complementing our previous work [SHT] [21 EH [33] the mathematical 
method presented here might be further investigated in several directions within the realm 
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of complexity science [III15H52]. In a previous work [7j we have introduced the concept of con¬ 
ditionally predictable dynamical systems which, in brief, can be stated as those systems whose 
trajectory can he exactly predicted provided that they met no exception in tracing it. In practice, 
a model for a conditional dynamical system is composed of a predictable part and a finite set 
of exceptions [7] that are introduced by means of suitable ;B-functions urn- Thus, the trajec¬ 
tory is explicitly known provided that no exception has been met on it. This, in turn, implies 
checking the whole trajectory for the occurrence of exceptions and, therefore, an associated 
degree of unpredictability and disbelief involved in the frailty of our prediction of the trajectory 
(if we renounce to carry out this check). By help of the B-function, the semipredictable dy¬ 
namical systems presented in this paper admit a straightforward generalization to conditionally 
semipredictable systems when exceptions are added as a function of the dynamical states reached 
by the trajectory on each point in the discrete spacetime. The layers can also be coupled in 
complex manners by adding suitable B-functions and realistic models for complex systems can 
be so designed for specific needs. 

• Music theory: CA have previously been used in algorithmic music composition |53[ [STj and we 
briefly sketch some ideas of how the methods presented in this article may contribute musical 
richness to this approach. In the theory of counterpoint [55] each of the different voices (instru¬ 
ments) evolves on a finite set of pitches forming an interesting melody that interweaves with the 
others. Each independent voice might be modeled by means of a suitable layer CA (probably 
of second or third order in time). Harmonic progressions, i.e. the sequences obtained from the 
evolution of simultaneous arrangements of the voices into chords, might thus be modeled by 
musically interesting graded CA rules. We conjecture that the latter might be constructed so 
as to evolve on Tymoczko orbifolds |56| thus warranting an efficient voice leading. A spatiotem- 
poral evolution of such a CA would thus yield in parallel an ensemble of musically consistent 
trajectories on these geometric structures. 

• Biophysics and medical physics: CA rules are often used as models for tumor growth and the 
spread of viruses and other diseases (see m and references therein). For example, in a very 
recent article m a 2D stochastic CA model with alphabet size p = A for the spread of the 
Ebola virus is considered. Although the CA there considered is stochastic, it is constructed out 
of three deterministic CAs working on four symbols and, hence, each being decomposable into 
two layers with po = pi = 2 symbols each. Decomposable CA rules might be useful not only 
to get insight into realistic dynamical models, but also to help in the design of specific drugs 
in order to target specific layers of a living organism afflicted by a malignant agent, by leaving 
other layers unaffected, and to elucidate the way in which the malignant agent (which itself 
might be modeled by a p-decomposable CA rule) affects the organism as a whole. 

• City landscapes and urban development: Cities might be seen as complex multilayered structures 
[581 159| emerging as a result of graded CA rules, some of them hierarchical. The mathematical 
methods presented in this paper, together with some others employed in the design of ‘ex¬ 
otic’ fractal surfaces IMIIMI Ea ED, might be applied to the understanding of these complex 
dynamical systems. 

We note that all results presented in this paper can be extended to two and three dimensions and 
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that this aspect does not at all affect the conclusions on the decomposability of the CA rules; All these 
higher dimensional CAs can be described by means of suitable extensions of Eq. by appropriately 
addressing the first argument of the digit function in Eq. for the problem at hand. We also note 
that layer CAs within a graded CA rule can also have different ranges (i.e. parameters I and r of the 
rules may be h dependent). The only parameter that is relevant to warrant decomposability is the 
size of the alphabet p. After appropriate straightforward modifications, the CAs to be decomposed 
can also be of higher order in time. 

We believe that the study of non-decomposable CA with p a prime number (to which all Boolean 
CA with p = 2 belong) is of fundamental interest, since they are the building blocks of all other CA 
dynamics in quite a similar way as prime numbers are building blocks for the non-negative integers. 
Further mathematical insights in non-decomposable CAs can be helpful to construct more subtle and 
increasingly realistic models of complex systems by means of graded CA rules, in particular, and 
p-decomposable rules, in general. 
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